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. Abstract. Quarkonia provide a sensitive probe of the properties of the hot dense 

| medium created in high energy heavy ion collisions. Hard scattering processes result 

£NJ ' in the production of heavy quark pairs that interact with the collision medium 

during hadronization. These in-medium interactions convey information about the 
fundamental properties of the medium itself and can be used to examine the 
■ modification of the QCD confining potential in the collision environment. Baseline 

C*~) \ measurements from p+p and d+Au collision systems are used to distinguish cold 

nuclear matter effects while measurements from heavy ion collision systems are used to 
, quantify in-medium effects. The PHENIX experiment has the capability of detecting 

' heavy quarkonia at 1.2 < |?7| < 2.2 via the [J, + fJ.~ decay channel and at |ry| < 0.35 via the 

e + e" decay channel. Recent runs have resulted in the collection of high statistics p+p 
data sets that provide an essential baseline reference for heavy ion measurements and 
allow for further critical evaluation of heavy quarkonia production mechanisms. The 
latest PHENIX results for the production of the J/ip in p+p collisions are presented 
^ ' and future prospects for ip' , \c and T measurements are discussed. 
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1. Introduction 

The study of quarkonia in p+p collisions provides a useful tool for probing heavy 
quark production mechanisms. Heavy quarks are predominantly generated in hadronic 
collisions via hard processes involving gluonic diagrams, but the details of hadronization 
remain unclear. Several theoretical models, with varying degrees of success, have been 
proposed to explain how a heavy quark pair evolves to form a bound quarkonium state. 
The color-singlet model generates J/ip particles from color singlet cc pairs that are in 
the same quantum state as the final J/ip PP- However, the J/ip production cross section 
is underestimated by an order of magnitude [2]. The NRQCD color-octet model [3] 
includes color-octet state cc pairs in addition to color-singlet state pairs that radiate soft 
gluons during hadronization to form a J/ip. Contrary to expectation, the color-octet 
matrix elements, derived from experimental data, are found to be non- universal [I]. 
Additionally, a large transverse polarization is predicted at high pt that is not observed 

| For the full list of PHENIX authors and acknowledgements, see Appendix 'Collaborations' of this 
volume 
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in the data. The color evaporation model provides a phenomenological approach to 
hadronization [5]. The charmonium states are formed in proportions determined by 
experiment for any cc pair below the DD threshold and hadronization occurs through the 
emission of soft gluons. To date, the most successful model of quarkonia production is 
a perturbative QCD approach involving 3-gluon mechanisms which is able to reproduce 
both the experimentally observed J/ip cross section and polarization [6]. Within 
this model, the hadronization mechanism includes a channel involving the fusion of 
a symmetric color-octet state with an additional gluon. The measurements reported in 
this paper provide further tests for these production mechanisms. 

2. PHENIX Experiment 

The PHENIX experiment is designed to detect heavy quarkonia at forward rapidity 
(1.2 < \r]\ < 2.2) via the decay channel and at mid-rapidity (\r)\ < 0.35) via the 

e + e~ decay channel jTJ. At forward rapidity the muon detectors are composed of cathode 
strip tracking chambers in a magnetic field and alternating layers of steel absorber and 
Iarocci tube planes that allow muons to be tracked and reconstructed over an acceptance 
of A0 = 360°. At mid-rapidity the central arm tracking detectors, composed of drift 
chambers, ring imaging Cerenkov detectors, and electromagnetic calorimeters, are used 
to detect electrons in two arms each covering A0 = 90° in azimuth. During the 2005 
run of the Relativistic Heavy Ion Collider the PHENIX experiment sampled 3.8 pb^ 1 of 
p+p collisions. From this data set a total of 8000 (1500) J/ip have been reconstructed 
in the dimuon (dielectron) channel [8]. 

3. Transverse Momentum Dependence of the J/tp Cross section 

As shown in Fig. (TJ the transverse momentum dependence of the J/ ip cross section can 
be mapped to 9 GeV/c using the new high statistics data. A comparison of the spectral 
shapes at mid and forward rapidity reveals that the pt spectrum is softer at forward 
rapidity. This is likely a result of the increased longitudinal momentum at forward 
rapidity, consequently there is less energy available in the transverse direction. 

The spectra are fit with the functional form, A x (1 + (pt/B) 2 )~ 6 , to extract the 
{pip). At mid-rapidity (p?) = 4.14 ± O.IS+q^o (GeV/c) 2 and the \ 2 P er degree of 
freedom (x 2 /ndf) is 23/19. At forward rapidity {p 2 T ) = 3.59 ± 0.06 ± 0.16 (GeV/c) 2 and 
the x 2 / n df is 28/17. If the exponent in the fit function is allowed be a free fit parameter 
a slightly better fit is achieved at forward rapidity, xV n df = 20/16, but the (p\) is 
not significantly modified, (p|>) = 3.68 (GeV/c) 2 . The previous results from the 2003 
run j9] yielded a significantly lower (p^) at forward rapidity, even though within errors 
the pt spectra agree. The increased statistics of the 2005 data allow for an improved 
understanding of the shape of the px spectrum due to the extended range in p T and 
the finer binning at low pt- The 2003 p+p results have been revisited and it was found 
that the systematic error was underestimated. A reanalysis of the 2003 run d+Au data 
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Figure 1. Left: The J/ip differential cross section times di-lepton branching ratio 
versus pt at mid and forward rapidity. The vertical error bars are the statistical 
and point-to-point uncorrelated error and the boxes are the point-to-point correlated 
systematic error. The solid lines are the fits [5]. Right: The J/ip (pi^) truncated at 5 
GeV/c plotted versus N part in p+p, Cu+Cu and Au+Au collisions [51 ITU1 ITT] . 

is underway to determine how this effects the interpretation of that data set. 

A comparison of the (pj) found in p+p collisions with that found in Au+Au and 
Cu+Cu collisions as a function of centrality is shown in Fig. [JJ [JJJl [8j [UJ. Because the 
heavy ion pr spectra only extend to 5 GeV/c the calculation of the {p\) is truncated at 
5 GeV/c for all data sets. Good agreement is found between the peripheral heavy ion 
data and the p+p measurements. 

4. Rapidity Dependence of the J/ip Cross section 

The rapidity dependence of the J/ip cross section has been mapped over the range 
—2.2 < y < 2.2, Fig. [2j The data exhibit a slight flattening over the rapidity range 
\y\ < 1.5. However, the systematic errors on the mid and forward rapidity points are 
independent, thus a narrower distribution is not excluded. A comparison of the shape 
of the distribution with various parton distribution functions and model calculations 
P, [121 [13] shows that many of the models, including the perturbative QCD approach 
involving 3-gluon mechanisms, are inconsistent with the steepness of the slope at forward 
rapidity and the slight flattening observed at mid-rapidity. 

5. Other Quarkonia 

In addition to the J/ip, the PHENIX detector has the capability of measuring the T, \ c 
and ip' quarkonia states in p+p collisions. Preliminary PHENIX measurements of the 
T as a function of rapidity and collision energy are available [Hj. The analysis of the 
Xc and ip' states is still ongoing. The 2006 p+p run provides a factor of 3 improvement 
in the statistics available for these measurements and the results will be forthcoming. 





Figure 2. J/ip cross section in p+p collisions as a function of rapidity [5] plotted with 
various parton distribution functions and model calculations for comparison [SI II 21 ITS'] . 



6. Summary 

With the existing high statistics data the PHENIX experiment has entered an era 
of precision J/ip measurements. The p+p data exhibit several features that provide 
challenges for production models. The pt spectrum has been measured over the range 
< px < 9 GeV/c and is softer at forward rapidity than at mid-rapidity. Furthermore, 
an improved assessment of the J/ip (p?) is possible with the 2005 data set and good 
agreement is found between the peripheral heavy ion data and the p+p measurements. 
The J /ip cross section as a function of rapidity is slightly flat and falls off rapidly 
at forward rapidity. This behavior is not well reproduced by the available model 
calculations. Future PHENIX data will shed light on these processes and open additional 
exciting avenues of quarkonia measurements. 



References 



[1] R. Baier et at, Phys. Lett. B102, 364 (1981). 
[2] F. Abe et al, Phys. Rev. Lett. 79, 572 (1997); 

[3] G. T. Bodwin et al, Phys. Rev. D51, 1125 (1995); erratum Phys. Rev. D55, 5853 (1997). 

[4] J. K. Mizukoshi, SLAC-PUB-8296, |hep-ph/9911384| (1999). 

[5] R. Gavai et al, Int. J. Mod. Phys. A10, 3043 (1995). 

[6] V. A. Khoze et al. , Eur. Phys. J. C39, 163 (2005). 

[7] K. Adcox et al, Nucl. Instr. Meth. A499, 469 (2003). 

[8] A. Adare et al, |hep-ex/0611020| (2006). 

[9] S. S. Adler et al., Phys. Rev. Lett. 96, 012304 (2006). 



10] A. Adare et al, nucl-ex/0611020 (2006). 



[11] H. Pereira Da Costa, Nucl. Phys. A774, 747 (2006). 

[12] R. Vogt, Heavy-Ion Phys. 25, 97 (2006) and private communication. 

[13] E. L. Bratkovskaya et al, Phys. Rev. C67, 054905 (2003). 

[14] M. Leitch, J. Phys. G32, S391 (2006). 



